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Abstract 

It has been conjectured by Ambrosanio, Kane, Kribs, Martin and Mrenna 
(AKM) that the CDF event pp — > e + e~77^T is due to a decay chain involving 
two neutralino states (the lightest and the next-to-lightest ones). The light- 
est neutralino (xakm) has been further considered by Kane and Wells as a 
candidate for cold dark matter. In this paper we examine the properties of 
relic xak m' s m their full parameter space, and examine the perspectives for 
detection by comparing theoretical predictions to sensitivities of various ex- 
perimental searches. We find that for most regions of the parameter space the 
detectability of a relic xakm would require quite substantial improvements 
in current experimental sensitivities. The measurements of neutrino fluxes 
from the center of the Earth and of an excess of p/p in cosmic rays are shown 
to offer some favorable perspectives for investigating a region of the xakm 
parameter space around the maximal tan (3 value allowed by the model. 
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I. INTRODUCTION 



The occurrence at CDF of the single event pp — > e + e~77pT [0] has prompted two different 
super symmetric interpretations 0-0, although also a non-supersymmetric explanation has 
been proposed 0. 

The two supersymmetric interpretations have a common scheme to explain the CDF 
event: first, an e + e~ pair is produced, pp — > e + e~, then the following decay chain takes 
place e ± — > e ± X 2 ,X 2 — > Xij, where X\ and X 2 are the lightest and the next-to-lightest 
supersymmetric particles (LSP and NLSP, respectively). The two supersymmetric inter- 
pretations differ in the identification of X\ and X 2 with definite supersymmetric particles. 
In one interpretation 0,0] X\ and X 2 are identified with the lightest and next-to-lightest 
neutralinos Xi an d X2, respectively; in the second interpretation |2|-|6| X 2 is the lightest 
neutralino Xi, whereas X\ is identified with the gravitino, which has the role of the lightest 
supersymmetric particle. 

To test whether or not one of the supersymmetric interpretations is correct, various pos- 
sible processes which can occur at CERN LEP or at Fermilab Tevatron have been discussed 
in Refs. 

Furthermore, should one of the supersymmetric interpretations be valid, this would have 
implications for the presence of supersymmetric particles as relics in the Universe; these 
could also provide a substantial contribution to the cosmological matter density 0. Thus the 
natural question arises, whether any of these fossil particles could be detected either directly 
or indirectly (or whether they are already excluded in force of the present experimental 
bounds). 

The case of relic gravitinos would practically represent a hopeless situation for experi- 
mental investigation, since gravitino interactions are too weak to allow detection jL(J. In 
models of gauge-mediated supersymmetry breaking (where the gravitino is the LSP) also 
the lightest messenger particle may be a viable dark matter candidate with a substantial 



relic abundance. This possibility has been recently investigated in Ref. |TT|. We do not 
pursue the discussion of the gravitino case any further here. 

Instead, in this paper we address the problem of the possible detection of relic neutralinos, 
whose specific properties are appropriate to a correct interpretation of the CDF event ||. 

As discussed in Refs. ||[7|], the interpretation of the CDF event in terms of the decay 
chain pp — > e + e~ , — > e ± X2, X2 XiT) se ^ s a number of very stringent constraints on the 
supersymmetric parameter space, and in particular on the nature of the neutralinos. These 
constraints, due to the kinematics of the event, and to the required sizes for the relevant 
cross section and decay branching ratios, imply that xi an d Xi are a very pure higgsino 
and a very pure photino, respectively. Detailed descriptions of the resulting supersymmetric 
parameter space are given in Ref. 0, and some of these results will also be reported here 
in the next section. For the moment, let us just anticipate some of the most prominent 
features of the model. The parameter tan /3 is in the very low side (i.e. 1 < tan f3 < 3) of its 
natural range: 1 < tan/? < 50; also, for the soft-breaking gaugino masses one has Mi ~ M 2 , 
rather than the usual relationship Mi ~ M2/2, motivated by unification assumption at 
Mqut (definitions and conventions for the supersymmetric parameters are as in Ref. |fL2|1 ). 
Furthermore, the kinematics of the CDF event (combined with the lower bound from LEP 
data) entails that m X2 — m Xl > 30 GeV and 30 GeV < m Xl < 65 GeV. In the following, a 
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neutralino eigenstate of the lightest mass, xi> wm be denoted by xakm, when its properties 
are those required by the supersymmetric interpretation of the CDF event, as suggested and 
described in Refs. [J|,[7j (LSP neutralino scenario). 

In ref. || some properties of a Xakm neutralino as a candidate for Cold Dark Matter 
(CDM) were discussed, and the perspectives for a direct detection were analysed, under 
the hypothesis that the contribution of Xakm to the cosmo logical density fl is substantial. 
The analysis was pursued there in the extreme case of a pure higgsino composition and in 
general for a parameter space sizably narrower than the one allowed by the supersymmetric 
interpretation of the CDF event. In particular, it was concluded that the perspectives for a 
direct detection for this candidate are rather favorable. 

In the present paper we reconsider the properties of Xakm as a relic particle, by ex- 
panding the previous analysis in many instances i) we explicitly take into account a 



possible gaugino-higgsino mixing in Xakm, which, although very tiny, might nevertheless 
have sizeable consequences in some processes for tan f} ~ 3; ii) we relax the requirement that 
Xakm contributes to Q significantly, since we wish to fully explore the experimental chances 
to detect a relic Xakm, even in the case it is not the main component of Cold Dark Matter 
(CDM); iii) we implement all the constraints from accelerators, including b — > s + 7 and the 
new bounds implied by the LEP2 measurements [TJ]]; iv) we examine what are the chances 



to detect relic Xaktm's using various detection strategies (direct detection as well as indirect 
measurements: neutrinos from the Sun |15 and from the Earth and the antiproton/proton 
ratio in cosmic rays). 

The motivations for the previous points are the following. On very general grounds 
one expects that a higgsino-like neutralino, such as Xakm, provides a large relic abundance 
fl x h 2 , but has very little chances to be detectable. Indeed, a higgsino interacts with matter 
through spin-dependent effects, whereas the sensitivities of the experimental searches which 
are based on neutralino-matter scattering, however expected to substantially improve in 
the near future, will still remain for a while only at the level of the much larger coherent 



effects |p!6| . Thus the perspectives for detection of a relic Xakm in the near future appear 
to be rather gloomy. Nevertheless, the conjecture of Refs. |3],[7||| is very challenging and 
potentially so much far-reaching, that it deserves a more careful analysis from the point of 
view of the actual perspectives of detectability. 

Therefore we have undertaken the present analysis, with the aim of investigating the 
various circumstances which could provide some better perspectives for experimental explo- 
ration of at least some physical region allowed to Xakm- This is why first, by taking into 
account the gaugino-higgsino mixing, even if small, we explore the possibility that coherent 
effects may help in providing direct and indirect detection rates with more substantial con- 
tributions than the ones due to the spin-dependent effects. Of course, this cannot occur for 
configurations with tan (3 ~ 1, but could happen for neutralino compositions at the upper 
side of the allowed range for tan /3, i.e. tan/3 ~ 3. Secondly, we have relaxed the constraint 
that Xakm provides substantial relic abundance, which in itself sounds rather arbitrary and 
at the same time forces Xakm to stay in a region of the parameter space where neutralino 
cross sections are automatically small. After all, should one be able to detect relic neutrali- 
nos compatible with the unique CDF event, this would already be a major breakthrough, 
even if these neutralinos do not provide a large Q\ Finally, apart from the more standard 
detection techniques for WIMPs (direct detection and detection of neutrinos from macro- 
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scopic bodies), also the antiproton/proton ratio in cosmic rays has been considered. Indeed, 
the Xakm neutralino holds some features which could favor this kind of signal: a small mass 
and a very tiny (but not vanishing) mixing [17 . 

A few more comments are in order here. All experimental data have to be correctly 
implemented in shaping the allowed parameter space for xakm- Therefore we have taken 
into account the b — > s + 7 process which is a very constraining bound to be implemented 
in any realistic model (it is not clear whether or not it was properly taken into account in 
the previous analyses of Refs. P,[15|l). Furthermore we have included in our analysis the 
very recent data from LEP2 . In the next section we show how these new results further 
constrain the parameter space of Refs. 0,0]. The bounds on the Higgs masses, which are 
important for our evaluation of the detection signals, are obtained from the experimental 
data of Refs. 

In section 2 we define the supersymmetric 
In sections 3 and 4 we 



The scheme of this paper is as follows, 
parameter space and examine some general properties of Xakm- 
present our results for the Xakm relic abundance and for its detection rates, respectively. 
Conclusions are finally reported in section 5. 



II. MODEL PARAMETER SPACE 

Our parameter space has been modelled according to the one of Ref. 0. It issues from the 
requirement that the CDF event is due to the process: pp — > e + e~, e^ 1 — > e ± X2,X2 — > Xi7- 
We only consider the case of ei production, which appears to be the favourite scheme among 
those suggested in Ref. M. 

The parameters are: M 1; M 2 , /z, tan/3, tua (mass of the CP-even Higgs neutral boson), 
mi = rriq (this is the common mass for sleptons and squarks, taken to be degenerate, with 
the exception of the left-handed selectron (of mass m e -J and of the lightest stop (of mass 
m^) and % (mixing angle in the stop mass matrix). 

Our analysis of the properties of a relic Xakm has been performed by varying the super- 
symmetric parameters of a low-energy Minimal Supersymmetric extension of the Standard 
Model (MSSM) in the following ranges |7j: 

region A 





1.05 


< tan/3 < 1.5 


55 


GeV 


< M 2 < 90 GeV 




0.8 


< M 2 /Mi< 1.2 


-70 


GeV 


<H < -33 GeV 


75 


GeV 


< m~ eL < 140 GeV 




m A 


60, 100, 200, 400 GeV 




rriq 


250, 500, 1000 GeV 


150 


GeV 


< m- tl < rriq 




-tt/2 


<Q~t <tt/2 



region B 
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1.5 < tan/? < 2.8 

40 GeV < M 2 < 130 GeV 
1.2 < M 2 /M 1 < 2 

-70 GeV < // < -33 GeV 

75 GeV < m iL < 140 GeV 

m A = 60, 100, 200, 400 GeV 

rriq = 250, 500, 1000 GeV 

150 GeV < m h < m- q 

-7r/2 <0 t - <vr/2 (2) 

In both cases mg l has been further required to satisfy the kinematical constraints among 
m e -j and the neutralino mass eigenvalues [JjJ. 

As usual, any neutralino mass-eigenstate is written as a linear superposition 

Xi = a il + ^ + c iH s + diH a (3) 

where 7, Z are the photino and zino states and H s , H a are defined by H s = sm(3H°+cos{3H%, 
H a = cos/9i?° — sin/5 if^? i n terms of the higgsino fields H%, supersymmetric partners 
of the Higgs fields H^, H%, which provide masses to the down-type and up-type quarks, 
respectively. 

The lightest neutralino state xii obtainable by varying the supersymmetric parameters 
in the regions A and B is what we define as a xakm neutralino. Its mass turns out to be 
confined in the range: 30 GeV < m x < 65 GeV (in region B the upper limit is about 60 
GeV). 

It is worth noticing here that in the previous analyses of relic Xakm's HIM], only a 
restricted region of the parameter space (tan/5 ~ 1) was considered, where higgsino purity 
in xakm is most pronounced. Also, in Refs. only the lowest part of the neutralino 

mass range was considered, 30 GeV < m x < 40 GeV, in order to avoid the Z-pole (and 
possibly Higgs-poles) in the neutralino pair-annihilation cross section, where the evaluation 
of the neutralino relic density requires great care. In this paper we include in our discussion 
both region A and region B of supersymmetric parameters. We also consider the whole 
m x range and discuss the effect on the detection signals of a careful calculation of the relic 
abundance over the poles of the neutralino-neutralino annihilation cross section. 

On the other side, as previously mentioned, the new constraints from LEP2 |14| have 



been included. We show some effects of these constraints in Fig.l. This figure displays in the 
plane fi — M2 those AKM configurations of regions A and B (see Eqs.(|l]||)) which survive 
the b — > s + 7 constraint. It turns out that some of them are already excluded by the new 
LEP data. The various curves denote the chargino isomass contours at fixed tan/5 which 



correspond to the current LEP lower bound on the chargino mass ||14|| . For a given value of 
tan p, the configurations on the right of the relevant line are disallowed. In particular, one 
sees that no AKM configuration survives for tan/5 > 2.6. 

The composition of Xakm is what establishes the size of the neutralino relic abundance 
and of the detection rates. This composition is shown in Fig. 2. In Sect. a of this figure we 
give the values of the weights |cti| 2 , |&i| 2 , |ci| 2 , |g?i| 2 of the 7, Z, H s , H a components for the 
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smallest value of tan (5: tan/3 = 1.05. In Sect.b we display the values of the same quantities 
cumulatively for the two values tan/3 = 2.15,2.5. We notice that, as anticipated, Xakm is 
largely dominated by H s , with a next-to-leading contribution from Z. In Sect.c of Fig. 2 
we give a scatter plot for the fractional weights of these two main components of Xakm 
over the full grid of Eqs.(|l], |2|). For some configurations of region B the value of the ratio 
|fri| 2 /|ci| 2 , even if small, may nevertheless be sizeable enough to allow coherent effects in 
neutralino-matter interaction to overcome the spin-dependent ones. 



III. RELIC ABUNDANCE OF xakm 

The neutralino relic abundance £l x h 2 is evaluated using the standard formula 

n x h 2 = 3.3 x IP' 38 , 1 (4) 



where 

r x f 

[Xi 



o 



dx < a ann v > . (5) 



< &ann v > is the thermally-averaged annihilation cross section times the relative velocity, 
g*{xf) is the number of degrees of freedom at the freeze-out temperature Tf and Xf — Tf/m x . 

Whenever the neutralino- neutralino annihilation cross section is not in the proximity 
of a pole or when anyway a great accuracy in the estimate of Q x h 2 is not important, we 
simply expand < a ann v > at small velocities < a ann v >= a + bx (x = T/m x ), and thus 
I(x f ) = ax f + bx f 2 /2 @J2T[ . 



Otherwise, when we are close to a pole for the annihilation cross section and we require a 
careful evaluation for the relic abundance, the function I(xf), which entails multiple integra- 
tions over x and over the two particle velocities, is carefully evaluated, in part analytically 
and in part numerically [^,^|. Since this procedure is much computer-time consuming, we 
have applied the following selection criteria. Out of the full set of neutralino configurations 
explored through a scanning of the regions A and B, we have selected a number of config- 
urations (denoted as set S in the following) which, according to our estimates of detection 
rates, have more chances to be detected in the future. The set S will be precisely defined 
later on. For the configurations of set S the relic abundance has been evaluated in the exact 
way (and compared to the approximated estimate), whereas for other configurations only 
the approximate method, based on the low-velocity expansion, has been adopted. 

In the evaluation of < cr ann v > all the // final states as well as the complete set of Born 
diagrams have been taken into account |[24f . 

In Fig. 3 we display the results of our evaluation. Diamonds and crosses represent the 
values of Q x h 2 for configurations of set S (diamonds denote the values of Q x h 2 calculated in 
the exact way, crosses give the values obtained with the low- velocity approximation). Dots 
denote the values Q x h 2 for the other configurations (in the low- velocity approximation). 

Some interesting features show up in this figure: i) Q x h 2 displays the typical dip at 
about 45 GeV (Z-pole); ii) in going through the pole in a ann , the approximated value of 
fl x h 2 changes from an overestimate to an underestimate of the correct value (see Sect.b of 
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Fig. 3); iii) as expected, in region A (tan/? ~ 1) the relic abundance may be quite sizeable, 
and may even fall in the favorite range QcDMh 2 — 0.2 ±0.1 [25], whereas in region B it 
turns out to be systematically below 0.03. 

The evaluation of Q x h 2 is important here not only to establish the role played by the 
Xakm neutralino as a CDM candidate, but also to provide the value of the local (solar 
neighborhood) density p x . This quantity enters in all the detection rates to be considered 
in the following. Here, to determine the value of p x , we adopt the following rescaling recipe 
26fl . For each point of the parameter space, we take into account the relevant value of 



the cosmological neutralino relic density. When Q x h 2 is larger than a minimal (Qh 2 ) min , 
compatible with observational data and with large-scale structure calculations, we simply 
put p x = pi. When fl x h 2 turns out to be less than (Qh 2 ) m i n , and then the neutralino may 
only provide a fractional contribution Q x h 2 /(Qh 2 ) m i n = £ to Qh 2 , we take p x = p£. The 
value to be assigned to (Qh 2 ) min is somewhat arbitrary, in the range 0.03 < (Qh 2 ) min < 0.2. 
In the present paper we have used (Qh 2 ) min = 0.03. As far as the value of pi is concerned, 
we have taken the representative value pi = 0.5 GeV • cm -3 . This corresponds to the central 
value of a recent determination of pi, based on a flattened dark matter distribution and 
microlensing data: pi = 0.5lto]i7 GeV • cm -3 p7 . 



IV. DETECTION RATES FOR XAKM 

The most natural question to be asked now is whether there may be some chance to 
detect a relic neutralino with the properties of Xakm- To provide an answer to this question 
we examine in detail three of the main methods for detecting relic particles (neutralinos in 
our case) [^8|: i) direct detection, ii) detection of neutrinos from macroscopic bodies (Earth 
and Sun), iii) measurement of an excess oip/p in cosmic rays, due to neutralino- neutralino 
annihilation in the halo. 

As mentioned in the introduction, on very general grounds one expects that the inter- 
action of the Xakm neutralino with matter takes place through spin-dependent effects 0. 
This is due to the fact that Xakm is an almost pure higgsino, and then couples to quarks 
mainly through a Z-exchange. This is particularly true for configurations where tan/3 ~ 1. 
However, as shown in Fig. 2, for xakm compositions at the upper extreme of the allowed 
tan/3 range, i.e. tan/3 = 2.5, the higgsino-gaugino mixing parameter |&i| 2 /|ci| 2 can reach 
a level of ~ 10% and then may switch on some coherent effects through Higgs-mediated 
or squark-mediated processes. This effect may trigger an enhancement in direct detection 
rates as compared to a simple evaluation based on spin-dependent effects only. A second 
beneficial effect due to a mixing in xakm is that also the neutrino outcome from the Earth 
(typically increased by coherent effects) might be sizably enhanced. 

A third detection method for xakm investigated in the present paper is the measurement 
of the antiproton component in cosmic rays. This experimental mean is very interesting, 
in view of the upcoming projects |29|j30[| , which should substantially increase the number of 
measured antiprotons, bringing the present total number of about 30 to something of order 
600 in a few- year time [3(J. This remarkable increase in statistics should soon allow us to 



discriminate between a fast- varying spectrum of secondary antiprotons and a flat spectrum 
of antiprotons of exotic origin for kinetic energies in the range 100 MeV < T < a few 
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GeV. Also from the theoretical point of view the peculiarity of xakm offers some interesting 
features for the p/p ratio (tiny higgsino-gaugino mixing and small mass). 

Now let us define our set S of AKM configurations. This is a set of representative points 
within regions A and B which satisfy the following prerequisites: their predicted signals 
either for detection of neutrinos from Earth and Sun (at least one of these) or for detection 
of antiprotons in cosmic rays is within two orders of magnitude from the current value of 
the relevant experimental upper bound (at 90 % C.L.). 



A. Direct detection 

Let us start our analysis of the detection methods from the most natural one: direct 
detection. This consists in the measurement of the energy released by a neutralino in its 
scattering off a nucleus in an appropriate detector, by using very different experimental 
techniques |2"B ]. Some of the most recent experimental results are given in Refs. |3"T|--|3"5]1 . 



In general, it is expected that, in the neutralino- nucleus scattering, coherent effects, when 
allowed by the neutralino composition, overcome spin-dependent effects. When this is the 
case, then the best way to compare experimental data, which usually refer to a variety of 
nuclear compositions, is to convert the upper limits on the energy spectra into upper bounds 
on the neutralino- nucleon scalar cross section (J^L ar - This procedure is a model-independent 
one, i.e. it does not depend on the neutralino composition. 

However, as previously discussed, for configurations of region A, the Xakm is a neu- 
tralino with a high higgsino-purity and then its spin-dependent interactions with matter are 
important. Therefore for configurations of region A the previous procedure is not the most 
appropriate one and consequently we consider a rate rather than the neutralino-nucleon cross 
section. As far as configurations of region B are concerned, it turns out that the maximal 
signals are already slightly dominated by coherent effects. Therefore, for our comparison 
between predictions and experimental upper limits for configurations of region B we use 
both quantities: rates and cross sections. 

For our evaluation of cross sections and rates for the process at hand, we used the method 



described in Refs. [p^j . Our results are reported in Figs. 4-5. In Fig. 4 we display the predicted 
values of the rate RNai for the scattering of a neutralino off a Nal detector, integrated over 
the range 3.75 KeV < E ee < 5.25 KeV, where E ee is the electron-equivalent energy. The 
reason for considering this quantity is that it provides one of the most stringent experimental 
upper bounds (for neutralinos interacting through spin-dependent effects and with a mass 



in the range 40 GeV < m x < 75 GeV): R^ai < 1 event/(Kg day) JHJ. From Fig. 4 it is clear 



that all the predicted values for AKM configurations fall far below the current experimental 
bound (by more than two orders of magnitude for both regions, region B being slightly 
better than region A). 

Fig. 5 displays the scatter plot of the neutralino-nucleon cross-section times the rescaling 
factor £ for configurations of region B only and compares these to the experimental upper 
bound. The experimental limit shown in this figure refers to an experiment using a Ge- 
detector | 35| . This bound is somewhat more restrictive than the previous one from the 



Nal-detector, since now we start dealing with coherent effects and then we can optimize 
all experimental data to obtain the most stringent upper limit. Consequently, the maximal 
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predicted signal for some configurations turns out to be a little closer to the current limit, 
but however away by about two orders of magnitude. We notice that some increase in R^ai 
and ^cr^liar * s due to the refined evaluation of the neutralino relic abundance (see Figs. 
4b-5). 



B. Neutrinos from the Earth and from the Sun 

Let us turn now to the possible signals consisting of fluxes of up-going muons through 
a neutrino telescope generated by neutrinos produced by pair annihilations of neutralinos 
captured and accumulated inside the Earth and the Sun. The evaluation of the muon 
fluxes, which is a rather elaborate multistep process, has been performed here according to 
the procedure described in Ref. [|57|] , to which we refer for details. 

In order to conform to the experimental data which we use as upper limits, we consider 
here fluxes of up-going muons integrated over muon energies above 1 GeV. The flux from 
the Earth ^ arth is also integrated over a cone of half aperture of 30° centered at the nadir, 
the one from the Sun $^ un is integrated over the whole possible outcome from the Sun, i.e. 
integrated over 25° around the Sun direction. We compare our evaluations to the Baksan 
upper limits: ^ arth < 2.1 x 10 _14 cni- 2 s- 1 (90%C.L.)> < 3 - 5 x 10~ 14 cm- 2 s- 1 (90%C.L.) 

Our results are shown in Figs. 6-7. We notice that for region A (Fig. 6a) the maximal value 
of $^" n , provided only by very few configurations at m x ~ 65 GeV, is ~ 5 x 10~ 15 cm _2 s _1 , 
anyway below the experimental upper bound roughly by a factor of 6. These configurations 
were disregarded in previous analyses [[],[nj. Most of the other configurations give signals 
largely spread over more than three decades. When we move from region A to region B 
(Fig. 6b), $^ un increases as expected, since coherent effects start playing some role in en- 
hancing the neutralino capture rate by the Sun. Here a significant number of configurations 
have a predicted level of $^ nra within an order of magnitude from the current experimental 
bound. For region B an even more favorable comparison between predictions and experi- 
mental sensitivity occurs for ^ arth (Fig. 7). Indeed, the maximal predicted value is away 
from the present experimental limit only by a factor 2. However, it is again apparent from 
Fig. 7 that the predicted values for ^ arth are spread over a very wide range of a few decades. 
As expected, the configurations with the highest values for the flux are those with a light 
Higgs boson A (m^ ~ 100 GeV) and some higgsino-zino mixing. The gap in between the 
two groups of configurations in Fig. 7 is indeed due to the step in rriA used in our sampling 
of the parameter space. 

We again notice that some increase in the level of the fluxes in Fig. 6b and Fig. 7 is 
due to the refined evaluation of the neutralino relic abundance. Furthermore it is worth 
noticing that some improvements in the comparison of the predicted values for ^ arth and 
the experimental data may be obtained through a more refined analysis of the fluxes in 
terms of their angular distribution |39| . 
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C. p/p in cosmic rays 



The annihilation of neutralinos in our halo may generate some amount of antiprotons in 
our Galaxy. A way of discriminating them, against the background due to the secondary 
antiprotons produced by primary cosmic particles with the interstellar medium, is to look 
at the p/p spectrum as a function of the kinetic energy T. At small T the p/p ratio due to 
secondaries increases quickly, as T increases, whereas the signal has a flat behaviour. 

Measurements of anti-protons have been going on for quite a long time with some con- 
flicting results [28|. More recent data 0,^31 seem to follow the behaviour of secondaries 
as evaluated in Ref. [52]. However, a much higher statistics is required to find out whether 
or not there might be a signal of some exotic origin for antiprotons. This very intriguing 
problem should be settled in a few-year time, due to upcoming experiments P5|J5D| which 



are expected to collect a total of about 600 antiprotons |3(J . 

We have evaluated the p/p ratio in the following way: i) the antiproton spectrum, as 
due to the neutralino pair annihilation in the halo, has been calculated as in Ref. |T7|], ii) 
the proton spectrum has been taken from Ref. [j43"|, iii) the propagation of the two fluxes 



has been evaluated using a leaky box model with an energy-dependent confinement time 
taken from Ref. 
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iv) the two spectra have been modulated by employing the procedure 
of Ref. [[EJ with the modulation parameter of Ref. [|T| and then integrated over the range 
250 MeV < T < 1000 MeV to conform to the experimental characteristics of one of the most 
significant experimental data: p/p = S.U±H x 10~ 5 for 250 MeV < T < 1000 MeV [jTJ. To 
make the comparison of our predicted values with the experimental data more meaningful, 
we use in the following the value p/p < 7.5 x 10" 5 as indicative of a 90 % C.L. limit for 
antiprotons of exotic origin. 

We show our results in Fig. 8. We notice that for a limited number of configurations 
in region B (Fig.8b) the predicted signal is rather close to the experimental value, but for 
many others the signals are away by orders of magnitude. The maximal predicted value 
for p/p is below the upper limit by a factor 3-4. Again it turns out that the improvement 
in the calculation of Q x h 2 enhances the expected signal. In Fig. 9 we give a scatter plot of 
p/p versus ^ arth to show how the same set of Xakm configurations provide the maximal 
predicted values for both of these two quantities. Fig. 10 shows where these configurations 
are located in the /i — M2 plane. 



V. CONCLUSIONS 

The conjecture that the CDF event pp — > e + e~ r yyfi T is due to a decay chain involving two 
neutralino states (the lightest and the next-to-lightest ones) ||[7| is certainly very intriguing, 
although great caution is in order, because of the existence of a single event of this sort and 
of the non-uniqueness in its interpretation Nevertheless, the interpretation of the 

event suggested in Refs. 0,[7|, if correct, would have so much impact on particle physics, 
that any possible experimental verification of it should be carefully investigated. Obviously, 
accelerators are the most suitable means for this purpose. 

Also the implications for relic supersymmetric neutralinos deserve much attention and 
experimental investigation. With this target in mind, we have extended the analyses of 
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Refs. p|,[15H in many ways. We have explored a much wider region in the supersymmetric 
parameter space than previously done and we have examined a variety of different detection 
means. In such a challenging enterprise of searching for a relic particle the only winning 
strategy is the one of combining as many independent searches as possible. 

Let us now summarize some of our results. We have considered three detection methods 
for relic neutralinos: i) direct detection, ii) detection of neutrinos from macroscopic bodies 
(Earth and Sun), iii) measurement of an excess of p/p in cosmic rays, due to neutralino- 
neutralino annihilation in the halo. For all of these search techniques we have evaluated the 
relevant signals exploring the widest allowed parameter space of the Xakm neutralino and 
we have compared our results to the current experimental bounds. In no case present ex- 
perimental data are such to provide information on a relic Xakm- All experimental methods 
require a substantial improvement in sensitivities before they may be capable of exploring 
some sizeable region of the Xakm parameter space. In general, the most easily accessible 
region is the one corresponding to the values of tan f3 close to the upper part of the AKM 
range 1.05< tan/? < 3 and to the smallest values of tua- Unfortunately, as discussed in 
Sect. 2, the new LEP2 data already exclude the Xakm configurations with tan/3 > 2.6. 

For direct detection the most easily accessible part of the Xakm parameter space requires 
a very significant experimental improvement in sensitivities of 2-3 orders of magnitude. 
We emphasize that our conclusion is based directly on consideration of neutralino-nucleon 
cross section, and then automatically takes into account the most stringent experimental 
measurements. 

The measurement of fluxes of upgoing muons from the center of the Earth appears to 
be in a much better situation, since some configurations are away from the present upper 
bound by a factor of two. However, an improvement of at least one order of magnitude in 
sensitivity would be necessary for an exploration of a significant number of configurations. 
The most suitable detector for this job appears to be MACRO with a muon energy threshold 
of about 1 GeV, whereas large-area neutrino telescopes such as AMANDA and NESTOR 
would not have much chances because of a much higher energy threshold. 

We have shown that some Xakm configurations may provide an excess of p/p in cosmic 
rays at a level that is away from the present measured value by a factor 3-4. This fact deserves 
much attention in view of the expected increase of statistics in the upcoming experiments 
in space. 

However, a word of caution is in order here. All our evaluations of signals for config- 
urations in region B are very sensitive to the value assigned to the parameter (Q x h 2 ) min , 
which enters in our rescaling of p x . Here we have used (fi x h 2 ) min = 0.03, which roughly 
corresponds to a minimal value for (fl x h 2 ) min . If a larger value of (Q x h 2 ) min = f x 0.03 
(/ greater than one) is employed, then in region B the estimates for R and $ M have to be 
reduced by a factor 1/ / and the values of p/p by a factor (l//) 2 . The same considerations 
apply to those configurations of region A, whose relic abundance is smaller than (f2 x /i 2 ) min . 

Finally, we point out that in the present paper our aim was to establish a comparison 
between the level of the maximal signals due to a Xakm neutralino and the experimental 
sensitivities currently available or expected in a near future. It is obvious that, in view of 
an actual experimental measurement, specific signatures and signal to background discrim- 
inations should be carefully investigated. This is beyond the scope of the present analysis. 
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Figure Captions 



Figure 1 - AKM configurations of regions A and B (see Eqs.([j],[D) which survive the 
b — > s + 7 constraint, displayed in the plane fi — M 2 . The various curves denote the chargino 
isomass contours which correspond to the current LEP lower bound on the chargino mass 
|14| , for different values of tan/3: tan/3 = 1.05 (solid), tan/3 = 1.5 (dotted), tan/3 = 2.6 
(dot-dashed), tan/3 = 2.8 (dashed). Dots (circles) denote configurations of region A (B). 

Figure 2 - Sect. a displays the values of the weights |ai| 2 , |&i| 2 , |ci| 2 , |<ii| 2 of the 
7,Z,i3 s ,i3 a components for the smallest value of tan/3: tan/3 = 1.05. Sect.b shows the 
values of the same quantities cumulatively for the values tan/3 = 2.15,2.5. Sect.c shows a 
scatter plot for the fractional weights of the two main components of Xakm over the full 
gridofEqs.(Hi). 

Figure 3 - Values of Q x h 2 for configurations in region A (B) are given in Sect. a (b). 
Diamonds and crosses represent the values of Q x h 2 for configurations of set S (diamonds 
denote the values of Q x h 2 calculated in the exact way, crosses give the values obtained with 
the low- velocity approximation). Dots denote the values Q x h 2 for the other configurations 
(in the low- velocity approximation only). The horizontal line corresponds to the value 
fi x h 2 = 0.03, below which we apply rescaling for the local neutralino density. 

Figure 4 - Predicted values of the rate RngI for the scattering of a neutralino off a Nal 
detector, integrated over the range 3.75 KeV < E ee < 5.25 KeV, where E ee is the electron- 
equivalent energy. Configurations in region A (B) are given in Sect. a (b). Diamonds and 
crosses denote the values of the displayed signal for configurations of set S when Q x h 2 is 
calculated with the exact expression and with the low-velocity approximation, respectively. 
Dots denote the values of the displayed signal for the other configurations, calculated in the 
low- velocity approximation only. 

Figure 5 - Scatter plot of the neutralino-nucleon scalar cross-section times the rescaling 
factor £ for configurations of region B. The line denotes the upper bound on the cross-section 
obtained using the experimental results of Ref. |}5| (Ge-detector). Diamonds, crosses and 
dots are as in Fig. 4. 

Figure 6 - The flux $^ un for configurations in region A (B) is given in Sect. a (b). 
Diamonds, crosses and dots are as in Fig. 4. The horizontal line denotes the Baksan upper 
limit: $J" n < 3.5 x 10- 14 cm- 2 s- 1 (90%C.L.) [§8J. 
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Figure 7 - The flux ^ arth for configurations in region B. Diamonds, crosses and 
dots are as in Fig. 4. The horizontal line denotes the Baksan upper limit: ^ arth < 
2.1 x 10- 14 cm- 2 s- 1 (90%C. J L.) g§. 

Figure 8 - p/p ratio in cosmic rays due to neutralino-neutralino annihilation in the 
galactic halo. The energy integration range is given in the text. The scatter plot refers to 
configurations of region A (B) in Sect, a (b). The horizontal line corresponds to the 90 % 
C.L. bound: p/p < 7.5 x 10~ 5 . Diamonds, crosses and dots are as in Fig.4. 

Figure 9 - p/p versus $^ art/l for configurations of region B. Diamonds, crosses and dots 
are as in Fig.4. The horizontal line corresponds to the 90 % C.L. bound: p/p < 7.5 x 10~ 5 , 
the vertical line denotes the Baksan upper limit: ^ arth < 2.1 x 10~ 14 cm~ 2 s~ 1 (90%C.L.). 

Figure 10 Location of the configurations of set S in the plane /i — M 2 . The curves 
are denoted as in Fig.l. Diamonds denote the configurations of set S, dots the other allowed 
configurations. 
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